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ABSTRACT 


An  analysis  of  an  experiment  to  determine  the  precision,  the  effects 
of  varying  amounts  of  environmental  data,  and  the  feasibility  of  the  use 
of  historical  environmental  data  at  the  NAVTORPSTA,  Keyport,  Washington, 
Dabob  Bay  facility  is  presented.   The  analysis  verifies  and  extends  the 
results  of  a  preliminary  investigation  of  the  experimental  data.   The  ex- 
periment consisted  of  tracking  a  75-kHz  source  at  six  depths  at  each  of 
six  horizontal  ranges.   The  data  were  analyzed  utilizing  an  isogradient 
ray-tracing  program.   The  results  showed  that  the  precision  of  the  range 
position  as  determined  by  the  ranging  system  is  dependent  only  on  the 
rise-time  of  the  75kHz  signal;   that  the  use  of  less  environmental  data 
and  historical  data  as  the  environmental  input  for  position  calculations 
have  been  shown  to  cause  only  a  minor  degradation  in  the  accuracy  of  the 
ranging  system. 
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I .   INTRODUCTION 

An  acoustic  experiment  was  conducted  to  provide  a  basis  for  the 
comparison  of  various  ray-tracing  techniques  which  can  be  utilized  to 
determine  the  three-dimensional  position  of  a  source  with  respect  to  a 
fixed  receiver.   The  experiment  was  conducted  at  the  Dabob  Bay  facility 
of  the  Naval  Torpedo  Station  (NAVTORPSTA) ,  Keyport,  Washington  on  21 
March  1974  by  members  of  the  Ray  Trace  Committee  of  the  Range  Study 
Group  at  the  Naval  Postgraduate  School,  Monterey,  California. 

The  purpose  of  the  experiment  was  to  determine  what  factors,  both 
instrumental  and  environmental,  influence  range  accuracy  and  precision. 
A  preliminary  investigation  of  two  of  the  fifty-seven  sets  of  data 
collected  [Ref.  1]  led  to  these  tentative  conclusions: 

(1)  The  precision  of  the  position  location  is  essentially  independent 
of  the  range  from  the  array  and  depends  only  on  the  rise-time  of  the 
75-kHz  signal. 

(2)  There  is  no  significant  difference  in  the  accuracy  of  results  ob- 
tained by  either  isovelocity  or  isogradient  ray-tracing  techniques. 

(3)  The  isogradient  method  appeared  to  need  less  environmental  data 
to  obtain  a  given  precision. 

This  thesis  was  undertaken  to  verify  and  to  extend  the  results  of 
this  preliminary  investigation  by  testing  the  isogradient  method  of 
ray-tracing  at  all  stations  and  depths  studied  and  by  determining  the 
amount  of  environmental  data  required  to  maintain  accuracy  and  by  ascer- 
taining if  historical  data  rather  than  real  time  data  may  be  utilized. 


II.   DABOB  BAY  TRACKING  RANGE 

A.   DABOB  BAY  FACILITY 

A  facility  for  the  three-dimensional  tracking  of  surface  ships, 
submarines,  and  torpedoes  is  located  at  the  Dabob  Bay  facility  of  the 
Naval  Torpedo  Station,  Keyport,  Washington.   Chosen  because  of  its 
favorable  oceanographic  features  and  its  proximity  to  the  NAVTORPSTA, 
Keyport,  Dabob  Bay  has  been  used  for  torpedo  testing  since  1949.   It 
consists  of  a  75-kHz  primary  tracking  system  enveloping  nearly  the 
entire  bay  providing  for  underwater  tracking  in  a  volume  approximately 
30,000  feet  by  4,500  feet  by  600  feet  deep. 

Dabob  Bay,  a  branch  of  the  Hood  Canal,  is  in  a  deep  depression  ad- 
jacent to  the  Olympic  Mountains.   The  Quilcene  River  flows  into  the 
northernmost  part  and  the  Dosewallips  River  into  the  southernmost.   The 
sides  of  the  bay  are  precipitous  and  predominately  rocky,  while  the 
bottom  is  mud.   Annual  precipitation  averages  about  15  inches  and  the 
average  snowfall  is  about  15  inches.   Winds  generally  blow  along  the 
length  of  Dsbob  Bay  in  a  Nortb-South  diiection.   Normal  wind  velocities 
are  5  to  15  miles  per  hour.   Tidal  levels  in  Dabob  Bay  range  from  about 
-5  to  +15  feet  referenced  to  mean  low-low  water  (MLLW) . 

Water  temperatures  change  markedly  from  season  to  season  and  even 
from  day  to  day.   Salinity  varies  considerably  with  depth  as  fresher 
surface  water  from  terrestrial  runoff  overlies  more  saline  water.   The 
variation  in  salinity  is  also  quite  seasonal,  reaching  a  maxiPium  in  the 
spring.   The  main  contributor  to  sound  velocity  changes  is  water  tem- 
perature variation  and  the  resulting  sound  velocity  profile  is 
representative  of  the  temperature  gradient. 


B.   RANGE  DESCRIPTION 

The  basic  underwater  tracking  system  is  comprised  of  three  compo- 
nents:  hydrophone  array,  computer  system,  and  vehicular  instrumentation. 
The  shipboard  transducer  emits  a  75-kHz  acoustic  pulse  in  synchronism 
with  a  master  clock  at  the  computer  site.   This  pulse  is  detected  by 
each  of  the  four  hydrophones  in  the  array  and  mixed  with  a  local  oscil- 
lator having  a  frequency  unique  to  that  hydrophone. 

The  outputs  of  these  oscillators  are  preamplified  and  fed  to  a 
multiplexer,  which  sums  the  signals  and  transmits  the  composite  signal 
via  underwater  cable  to  the  computer  site,  where  it  is  processed  to  de- 
termine the  in-water  transit   time  to  each  hydrophone  (Tx,  Ty ,  Tz,  Tc) . 

1.   Hydrophone  Array  Configuration 

The  hydrophones  in  the  NAVTORPSTA  Range  are  arranged  in  a 
short-base  line  system  in  groups  of  four,  each  one  located  on  one  of 
four  adjacent  corners  in  an  imaginary  30-foot  cube,  thereby  defining 
the  orthogonal  coordinate  system  in  which  the  measurements  were  made. 
These  hydrophone  arrays  use  buoyant  spheres  which  exert  an  upward  force 
in  excess  of  two  thousand  pounds  to  keep  them  upright.   The  array  is 
not  free  to  rotate  about  the  Z-axis ,  but  may  tilt  as  much  as  two  degrees 
from  the  horizontal  in  strong  currents.   This  tilt  is  measured  by  servo- 
pendulum  transducers  whose  specified  resolution  is  0.00833  degrees, 
accuracy  0.025  degrees,  and  linearity  0.05  per  cent  over  the  full  scale. 

The  range  consists  of  six  arrays  (Fig.  1)  aligned  along  a  line 
450  feet  east  and  parallel  to  the  range  axis,  spaced  at  6000-foot  in- 
tervals, with  the  00  array  being  the  farthest  north.   The  range  axis  is 
191°  18'  14"  T. 
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2.   Computer  System 

The  computer  system,  a  Scientific  Data  Systems  (SDS)  Model  920, 
is  installed  on  Zelatched  Point  and  consists  of  three  main  subsystems: 
the  signal-processing  subsystem,  the  data-collection  subsystem,  and  the 
computer  subsystem.   The  signal-processing  subsystem,  the  link  between 
the  hydrophones  and  the  data-collection  subsystem,  processes  the  multi- 
plexed signal  received  from  the  arrays,  discriminates  against  unwanted 
signals,  and  determines  the  in-water  transit  time  to  each  hydrophone. 
The  data-collection  subsystem  establishes  the  master  clock  timing  by 
which  all  timing  is  computed.   It  also  prepares  all  the  array  signals 
for  the  computer  subsystem,  which  calculates  the  position  of  the  tracked 
object,  prints  the  tracking  data,  and  records  the  data  on  magnetic  tape. 
For  the  real-time  examination,  the  data  were  plotted  on  X-Y  plotters. 
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111 '       EXPERIMENTAL  PROCEDURE 

A  92.5-dB  (re  1 /'bar)  acoustic  source  with  approximately  ten  pounds 
of  ballast  was  suspended  from  the  stern  of  a  sound  boat  to  depths  of  25, 
50,  75,  100,  150,  and  200  feet  at  six  horizontal  ranges  (stations)  from 
900  feet  to  4700  feet  from  the  02  hydrophone  (Fig.  2).   The  first  sta- 
tion was  located  approximately  900  feet  north  of  the  array,  and  250  feet 
east  of  the  range  axis.   Successive  stations  were  located  parallel  to 
the  range  axis;   the  second  being  directly  over  the  array,  the  third  600 
feet  past  it,  and  then  at  about  1500-foot  intervals  to  about  4500  feet, 
beyond  which  signal-to-noise  problems  made  further  examination  impossi- 
ble.  The  source  emitted  one  1.3  millisecond  75  kHz  pulse  per  second  for 
sixty  to  ninety  seconds  at  each  depth  at  each  station. 

The  horizontal  position  of  the  boat  on  the  range  was  found  from  the 
range's  Autotape  system,  a  portable  commercial  microwave  system  used 
exclusively  for  ship  tracking  that  yields  accuracies  of  +  3  feet.   The 
position  of  the  acoustic  source  with  respect  to  the  Autotape  antenna  was 
determined  from  the  length  of  cable  let  out,  the  horizontal  distance 
from  the  antenna  to  the  stern  of  the  boat,  and  the  ship's  heading. 
Tidal  data  were  recorded  and  used  for  correcting  depth  coordinates  to 
MLLW.   The  array  is  known  to  be  585  feet  below  MLLW,  so  that  the  depth 
of  the  source  below  the  surface  and  the  tide  measurement  allows  calcula- 
tion of  the  vertical  distance  between  the  source  and  the  array.   All 
measurements  were  made  with  respect  to  the  array's  coordinate  system, 
although  the  array's  orientation  and  location  on  the  range  were  known, 
making  it  possible  to  transform  this  data  to  a  true  position  on  the 
range.   The  X-axis  of  the  02  array  was  25°  24'  36.0"  T. 
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A  NAVTORPSTA  digital  velocimeter  that  has  an  accuracy  to  +1.0  ft/sec 
and  resolution  of  +  0.7  ft/sec  was  used  to  measure  the  sound-velocity 
profile  at  each  station.   NAVTORPSTA  instrumentation  at  Dabob  Bay  meas- 
ured and  recorded  the  acoustic  travel  times  from  the  source  to  each  of 
the  four  hydrophones  in  the  array  to  six  decimal  places  for  every  source 
depth  at  each  station.   The  three-dimensional  position  of  the  source 
with  respect  to  the  array  was  calculated  by  the  NAVTORPSTA  NUTRAK  III 
computer  program. 

The  data  were  gathered  on  21  March  1974  during  a  neap  tidal  period 
and  while  the  bay  surface  was  smooth  and  glassy,  although  there  had  been 
ripples  earlier  in  the  day.   The  day  was  bright  and  sunny  with  scattered, 
thin,  high  overcast,  the  temperature  was  55° F,  and  no  noticeable  wind  was 
present.   The  shiny  stainless  steel  transducer  disappeared  from  view  when 
about  eight  feet  deep,  indicative  of  thick  biological  concentrations  in 
the  water.   Rudimentary  drift  measurements  of  boat  movement  relative  to 
the  surface  water  were  taken  by  throwing  paper  into  the  water  and  watch- 
ing its  relative  motion  with  respect  to  the  boat,  which  indicated  a  value 
of  approximately  0.3  ft/sec. 
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IV .   RANGE  DETERMINATION  TECHNIQUES 

A.  TECHNIQUES 

The  three-dimensional  position  of  an  acoustic  source  with  respect  to 
a  receiver  may  be  readily  determined  by  computer  analysis.   Two  ray- 
tracing  techniques  which  are  easily  adaptable  for  high-speed  digital 
computation  are:   (1)  an  isovelocity  technique  in  which  the  water  mass  is 
divided  into  layers  of  isovelocity  water;   and,  (2)  an  isogradient  tech- 
nique in  which  the  water  is  divided  into  isogradient  layers.   Both 
techniques  utilize  Snell's  law,  which  describes  the  refraction  of  sound 
rays  in  a  medium  of  variable  velocity. 

The  isovelocity  method  is  currently  utilized  in  the  NAVTORPSTA  pro- 
cedure NUTRAK  III.   STUTRACK  I  was  developed  in  the  preliminary 
experimental  analysis.   It  consists  of  both  an  isovelocity  and  an  iso- 
gradient method  to  be  utilized  as  comparators  with  NUTRAK  III.   Both 
programs  require  as  inputs  the  velocity  profile,  the  transit  time  and 
angle  of  arrival  of  the  acoustic  wave  at  the  geometrical  center  of  the 
array. 

B.  PREPARATION  OF  DATA 

The  required  inputs  for  the  computer  analysis  must  be  derived  from 
the  information  available  at  the  array,  which  is  the  transit  times  of  the 
acoustic  pulse  from  the  tracked  source  to  the  four  hydrophones  of  the 
receiving  array.   The  geometry  involved  in  the  calculation  of  the  raw 
array  coordinates  is  shown  in  Figure  3.   The  expressions  for  the  raw 
array  coordinates  are: 
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RAWX  =  _^L   (Tc  +  Tx)   (Tc  -  Tx)  (1) 

2d 

2 

RAWY  =  -£t~   (T  +  T  )   (T   -  T  )  (2) 

2d     c    y     c    y 

and 

2 

RAWZ  =  -~-   (T  +  T  )   (T   -  T  )  (3) 

2d     c    z     c    z 

where 

c  =  the  average  sound  velocity,  ft/sec; 

d  =  the  distance  between  hydrophones,  30  ft; 

T  ,  T  ,  T  ,  T  =  travel  times  to  respective  hydrophones,  sec. 

A  correction  factor  must  be  applied  to  correct  for  the  array  receiver 
integration  time  (pulse  width  discrimination),  800  /is .   The  relations  now 

become: 

*-* 

(T«  +  T'  -  K)   (T'  -  Tf)  (4) 

ex         ex 

(T1  +  T1  -  K)   (T1  -  T')  (5) 

c    y         c    y 


and 


where 


0 

*-* 

RAWX 

— 

C 

2d 

RAWY 

= 

c2 

2d 
2 

RAWZ 

= 

c 

(T'  +  T*  -  K)   (T'  -  T')  (6) 

2d      c     z  c     z 


-3 

K  =  twice  the  integration  time,  1.6  x  10  s. 

T',  T1,  Tf,  T'  =  the  measured  transit  time,  sec. 
x   y   z   c 

These  raw  coordinates  are  referenced  to  a  non-horizontal  reference 
plane  due  to  the  array's  freedom  to  tilt  in  the  horizontal.   The  array  is 
instrumented  to  sense  the  angles  the  array  axis  makes  with  the  horizontal, 
XTILT  and  YTILT. 

The  correction  of  the  RAW  coordinates  is  a  matter  of  geometry,  as 
shown  in  Figure  4.   The  corrected  coordinates  become: 


15 


CORX  =  RAWX  -  RAWZ  [SIN(XTILT)] 
CORY  =  RAWY  -  RAWZ  [SIN(YTILT)] 
CORZ  =  RAWZ  +  RAWX  [SIN(XTTLT) ]  +  RAWY  [SIN(YTILT)] 


(7) 
(8) 
(9) 


In  a  homogeneous  bodv  of  water,  the  minimum  time  path  between  an  output 
transducer  and  a  receiving  hydrophone  is  a  straight  line.   When  the  speed 
of  propagation  is  spatially  varying  due  to  temperature  and  salinity  vari- 
ations, refraction  occurs  and  the  minimum  time  path  is  no  longer  a  straight 
line  (Figure  5) .   Corrections  must  therefore  be  made  to  determine  the  true 
position  of  the  source. 

From  Snell's  law,  a  minimum  time  path  is  computed,  working  back  in 
time  from  the  array  center  in  the  direction  of  the  "apparent"  position  and 
tracing  the  acoustic  ray  through  each  velocity  layer  until  the  measured 
time  is  consumed.   The  vertical  angle  of  entry  at  the  array,  computed  from 
the  tilt-corrected  array  coordinates  is 


A  =  SIN 
o 


-1 


CORZ' 


CORX2  +  CORY2  +  CORZ2 


1/2 


(10) 


This  entry  angle  is  assumed  to  be  equal  to  the  entry  angle  for  the 
minimum  time  path.   This  assumption  is  good  if  the  source  is  far  enough 
so  that  the  wave  fronts  are  planar  over  the  dimensions  of  the  array,  and 
the  speed  of  sound  is  constant  over  the  dimensions  of  the  array. 

The  travel  time,  T,  to  the  origin  is  computed  from 


(ID 


T       R 
c       c 


where 


S  -  [(CORX)2  +  (CORY)2  +  (CORZ)2]   1/2 


(12) 


and 
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R  -  [(CORX  +  15)2  +  (CORY  +  15)2  +  (CORZ  +  15)2]   1  (13) 

c 

where  R  is  the  slant  range  to  the  array  center  hydrophone,  and  S  is  the 

slant  range  to  the  array  center. 

The  entry  angle,  A  ,  the  travel  time,  T,  along  with  velocimeter  data 

describing  the  sound  velocity  profile  may  now  be  utilized  as  inputs  to 

the  computer  program.   The  outputs  will  provide  the  three-dimensional 

position  of  the  source  with  respect  to  the  receiver. 
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V .   PRELIMINARY  ANALYSIS 

A.  ANALYSIS 

The  analysis  of  the  experimental  data  was  to  have  determined  the 
accuracy  of  NUTRAK  III  as  compared  with  both  modes   of  STUTRACK  I.   The 
Autotape  system  was  to  form  the  reference  for  this  determination,  but  its 
precision  was  much  less  than  expected.   The  analysis  was,  therefore,  pri- 
marily devoted  to  the  determination  of  the  precision  of  the  various 
methods  of  ray- tracing.   Stations  1-25  and  6-50  were  chosen  for  prelimi- 
nary investigation  since  station  1  was  closest  to  the  array  and  station 
6  was  farthest.   The  effects  of  varying  the  amount  of  environmental  data 
was  also  investigated. 

B.  RESULTS 

NUTRAK  III  showed  good  point  to  point  continuity  at  station  1  for  all 
coordinates,  while  the  coordinates  at  station  6  were  very  erratic.   A 
transformation  of  the  coordinates  to  a  spherical  coordinate  system  showed 
that  the  variations  of  phi  versus  theta  to  be  approximately  the  same  at 
both  stations. 

Both  modes  of  STUTRACK  I  were  compared  to  NUTRAK  III  and  in  both 
cases  no  significant  differences  were  noted  when  the  most  structured 
sound  velocity  profile  was  utilized.   As  the  number  of  linear  segments 
defining  the  sound  velocity  profile  was  decreased,  the  isovelocity  tech- 
nique provided  consistently  poorer  results.   Using  the  isogradient 
technique  to  compute  the  source's  position  resulted  in  only  a  slight 
degradation  as  the  number  of  segments  in  the  profile  was  decreased  and  was 
considered  adequate  in  precision. 
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C.   CONCLUSIONS 

Comparison  of  the  two  stations  investigated  showed  station  6-50  to 
have  greater  spatial  variability.   The  transformation  to  spherical  co- 
ordinates shows  the  difference  is  not  present  in  the  angular  display  as 
would  be  expected  if  the  time  errors  were  the  important  source  of  vari- 
ability.  Improved  range  accuracy  can  only  be  obtained  by  improving  the 
accuracy  of  the  time  measurements  or  increasing  the  baseline  of  the 
hydrophone  array.   The  isogradient  method  of  STUTRACK  I  provides  greater 
precision  than  the  two  isovelocity  techniques  when  less  environmental 
data  was  utilized. 
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VI.   PRECISION  OF  DATA 

The  data  collected  at  Dabob  Bay  were  tested  at  all  stations  to  verify 
that  the  results  of  the  preliminary  analysis  were  generally  applicable. 
If  the  precision  was  independent  of  range  and  depth  the  standard  devia- 
tion of  the  source  coordinates  as  calculated  from  the  transit  times  at 
the  array  should  be  zero  with  a  small  standard  deviation  about  the  mean. 
Factors  which  could  cause  variation  are  the  rise-time  of  the  75-kKz  source 
(as  discussed  in  the  preliminary  analysis) ,  and  the  actual  drift  of  the 
source  transducer  caused  by  the  currents  which  exist  in  the  experimental 
area. 

Analysis  of  the  precision  of  the  ranging  system  must  remove  the  ef- 
fects of  drift  from  the  calculations.   Regression  analysis  on  the 
components  of  the  calculated  horizontal  position  of  the  source  provided 
an  equation  for  the  displacement  of  the  transducer  as  a  function  of  time. 
A  second-order  regression  equation  was  considered  to  adequately  describe 
this  displacement  as  it  accounts  for  the  velocity  and  acceleration  com- 
ponents of  the  drift. 

Each  set  of  data  was  processed  to  determine  the  calculated  components 
of  the  source  position.   Each  horizontal  component  was  fitted  to  a  simple 
second-order  curvilinear  regression  equation  with  time  being  the  indepen- 
dent variable  and  the  corresponding  component  as  the  dependent  variable. 
The  value  of  the  component  evaluated  from  the  regression  equation  for  each 
time  was  taken  as  the  estimated  value  at  that  time.   By  taking  the  dif- 
ference between  the  estimated  and  the  corresponding  calculated  values  the 
effects  of  drift  are  removed  from  the  computations.   The  mean  depth  was 
taken  to  be  the  simple  arithmetic  mean.   The  standard  deviation,  standard 
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error  of  estimation,  of  these  differences  between  the  estimated  and  cal- 
culated coordinates  is  a  measure  of  the  precision  of  the  ranging  system. 
Tables  I-III  show  the  standard  errors  for  the  theta,  phi,  and  radius 
components  of  the  source  position  for  all  stations  and  depths.   Figures 
6  and  7  are  scatter  diagrams  of  the  standard  errors  of  phi  versus  theta 
by  station  in  Figure  6  and  by  depth  in  Figure  7.   Both  Figures  show  that 
the  standard  errors  with  the  drift  effects  removed  are  essentially  random 
with  no  dependence  on  station  or  depth. 
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VII.   ENVIRONMENTAL  EFFECTS  ON  RANGE  PRECISION 

The  effect  on  range  position  accuracy  and  precision  due  to  varying 
the  amount  of  environmental  information  was  investigated  by  varying  the 
number  of  linear  segments  which  were  used  to  describe  the  sound  velocity 
profile.   Decreasing  the  number  of  segments  which  describe  a  particular 
profile  while  maintaining  the  required  precision  would  result  in  faster 
processing  time  for  the  ranging  data. 

Three  methods  of  decreasing  the  number  of  segments  in  a  profile  were 
studied:   (1)  equal  changes  in  gradient;   (2)  equal  changes  in  velocity; 
and  (3)  "eyeball"  fitting.   All  three  methods  were  compared  to  determine 
which  would  suffer  the  least  degradation  in  accuracy  and  precision  as 
the  profile  was  simplified.   In  the  experiment,  environmental  data  was 
provided  by  a  velocimeter  which  described  the  existing  sound  velocity 
profile  in  64  to  67  discrete  points  in  600  feet.   Any  profile  described 
adequately  with  fewer  points  would  significantly  reduce  data  processing 
time.   The  equal-change-in-gradient  and  velocity  methods  utilized  in  this 
investigation  were  developed  by  the  author  for  this  analysis. 

The  equal-change-in-gradient  method  determines  a  set  of  linear  seg- 
ments by  locating  the  depth  and  velocity  points  which,  if  connected  by 
straight  lines,  would  cause  a  given  percentage  change  in  the  gradient. 
The  first  three  points  from  the  bottom  of  the  most  detailed  sound  velo- 
city profile  are  fitted  with  a  straight  line  by  the  method  of  least 
squares.   The  slope  of  this  line  is  used  as  the  initial  reference  gradient 
The  gradient  of  the  linear  segment  in  the  original  profile  between  the 
first  point  used  to  determine  the  reference  gradient  and  the  next  point 
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in  the  profile  is  computed.   This  gradient  is  compared  with  a  test  value 
given  by  the  reference  gradient  plus  a  designated  change  from  5  to  30%. 

When  the  computed  gradient  exceeds  the  test  value,  the  prior  point  in 
the  profile  is  considered  as  the  end  point  in  that  segment.   The  points 
thereby  included  in  this  segment  are  fitted  by  the  method  of  least  squares. 
The  end  point  of  the  segment  which  has  been  processed  is  now  the  first 
point  used  in  the  determination  of  a  new  reference  gradient  as  described 
above  and  the  next  linear  segment  is  determined  by  repeating  the  process. 
The  velocity  at  the  highest  depth  of  the  initial  segment  does  not  equal 
the  velocity  at  the  lowest  point  of  the  second  segment  after  processing; 
therefore,  an  adjustment  must  be  made  to  cause  the  two  segments  to  inter- 
cept at  this  mutual  depth.   This  adjustment  consists  of  making  the  sound 
velocity  of  both  profiles  the  mid-point  of  the  calculated  velocities  at 
that  depth.   This  process  is  continued  until  the  surface  is  reached. 
Appendix  A  is  a  listing  of  the  computed  program  which  utilizes  this 
method  for  decreasing  the  number  of  segments  which  describes  a  sound 
velocity  profile. 

The  equal-change-in-velocity  method  utilizes  a  constant  percentage 
change  in  velocity  between  the  minimum  and  maximum  velocities  in  the  pro- 
file and  the  minimum  and  maximum  velocities  to  decrease  the  number  of 
segments  in  the  profile.   The  minimum  and  maximum  sound  velocities  of  the 
original  profile  are  located;   a  fixed  percentage  (5  to  30%)  of  the  dif- 
ference between  them  is  used  to  determine  the  velocity  at  each  point  in 
the  new  profile.   The  surface  sound  velocity  is  incremented  by  the  given 
change  in  velocity  and  the  original  profile  is  searched  until  the  depth 
of  the  incremented  velocity  is  located  or  the  minimum  or  maximum  velo- 
cities occur.   The  method  determines  the  direction  of  the  change  in 
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velocity  from  the  original  profile.   Once  the  new  point  in  the  profile  is 
determined,  its  velocity  is  incremented  and  the  sequence,  as  described 
above,  continues  until  the  bottom  of  the  original  profile  is  reached. 
The  original  sound  velocity  profile  is  now  described  by  a  lesser  number 
of  linear  segments.   Appendix  B  is  a  listing  of  the  computer  program 
which  utilizes  this  method. 

The  "eyeball"  method  is  simply  the  author's  estimate  of  the  best- 
fitted  segments  which  describe  a  sound  velocity  profile.   It  is  derived 
by  estimating  the  segments  which  provide  the  best  description  of  the 
original  profile  in  the  minimum  number  of  segments.   Then  by  decreasing 
the  number  of  possible  segments,  estimates  are  made  of  the  best-fitted 
profile  with  the  decreased  number  of  segments. 

Thirteen  curves  describing  the  sound  velocity  profiles  at  both  sta- 
tions 1  and  6  were  determined  using  the  above  methods,  six  by  the  equal- 
change-in-velocity  method,  three  by  the  equal-change-in-gradient  method, 
and  four  by  the  "eyeball"  method.   Figures  8-10  show  the  change  in  the 
characteristics  of  the  profile  at  station  1  due  to  its  description  by 
each  of  these  methods.   These  profiles  were  used  to  study  the  effect  on 
source  position  due  to  varying  the  amount  of  environmental  data.   The 
position  of  the  source  at  each  depth  for  stations  1  and  6  was  calculated 
for  each  profile.   Tables  IV-VII  show  the  resultant  standard  errors  of 
estimate  from  these  calculations.   At  both  stations  the  standard  error 
of  the  theta  and  phi  components  remained  constant  as  the  amount  of  envi- 
ronmental data  was  varied.   The  standard  error  of  the  slant  range  changed 
with  each  profile.   Excluding  the  results  for  the  "eyeball"  method  at 
station  6,  the  average  standard  error  for  all  methods  was  0.181  feet  at 
station  1  and  0.317  feet  at  station  6.   The  average  standard  errors  by 
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method  for  stations  1  and  6  were:   (1)  equal-change-in-velocity ,  0.160 
and  0.343  feet;   (2)  equal-change-in-gradient ,  0.220  and  0.318  feet;   and 
(3)  "eyeball",  0.172  and  1.586  feet.   The  equal-change-in-velocity  method 
appears  to  cause  the  least  degradation  in  positional  precision  with  the 
least  environmental  data  at  these  two  stations. 
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VIII.   HISTORICAL  ENVIRONMENTAL  ANALYSIS 

The  possibility  of  using  historical  environmental  information  rather 
than  real-time  information  as  the  sound  velocity  input  for  computing 
position  has  as  its  major  advantage  the  lessening  of  the  cost  and  time 
required  prior  to  performing  any  range  operation.   Historical  information 
along  with  its  reduction  by  the  methods  discussed  in  Chapter  VII  would 
decrease  total  processing  time  and,  thereby,  increase  the  efficiency  of 
the  ranging  system.   This  hypothesis  is  only  valid  if  the  use  of  this  in- 
formation does  not  degrade  the  accuracy  of  the  positional  calculations. 
The  objective  of  this  section,  therefore,  is  to  analyze  the  feasibility 
of  utilizing  historical  environmental  information  in  computing  position 
with  precision  as  the  measure  of  effectiveness. 

Figures  11  thru  22  show  the  twelve  sound  velocity  profiles  which  de- 
scribe the  environmental  conditions  which  exist  at  Dabob  Bay.   The 
profiles  were  based  on  three  years  (1971-1973)  of  environmental  informa- 
tion collected  at  the  facility  [Ref.  2].   Each  profile  type  shows  the 
minimum,  maximum,  and  median  (mid)  profiles  which  are  characteristic  of 
the  particular  type.   These  historical  sound  velocity  profiles  were  used 
in  this  analysis. 

In  order  to  perform  the  analysis,  the  inputs  for  the  range  positioning 
algorithm,  STUTRACK  I,  had  to  be  determined  based  on  the  historical  pro- 
files.  A  simulation  was  developed  to  generate  the  transit  time  of  an 
acoustic  pulse   from  a  source  to  the  array,   the  angle  of  entry  of  the 
acoustic  wave  with  respect  to  the  center  of  the  array,  and  the  apparent 
horizontal  distance  of  the  source  from  the  array.    The  simulation 
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generates  these  parameters  by  searching  for  the  entry  angle  in  which  a 
ray  would  enter  the  array  from  a  known  position  of  a  source.  Once  the 
entry  angle  is  located,  the  transit  time  for  that  angle  is  determined. 
This  simulation,  as  with  STUTRACK  I,  utilizes  Snell's  law  as  its  basis 
and  performs  the  same  calculations  as  STUTRACK  I,  but  in  reverse  order. 
Appendix  C  is  the  computer  listing  of  the  subroutine,  ANGCAL,  developed 
for  the  simulation. 

The  simulation  computes  the  input  data  for  STUTRACK  I  at  a  fixed 
horizontal  position  from  the  array.   This  position  was  chosen  to  provide 
a  horizontal  range  from  the  array  similar  to  that  of  station  6  in  the 
experiment.   Fifty  vertical  points  at  four-foot  intervals  from  the  sur- 
face to  200  feet  were,  used  as  the  vertical  component  of  the  source 
position.   The  use  of  this  scenario  was  chosen  to  provide  data  for  the 
analysis  under  the  conditions  that  were  found  to  cause  the  least  preci- 
sion of  source  position  in  the  experiment,  station  6. 

Each  of  the  historical  profiles  was   first  described  in  60  ten-foot 
layers  from  the  surface  to  600  feet,  and  this  description  was  used  as  the 
environmental  input  for  the  simulation.   The  transit  time,  entry  angle, 
and  apparent  horizontal  range  for  the  50  positions  were  used  to  compute 
the  three-dimensional  position  of  the  source  with  STUTRACK  I.   Each 
historical  profile  was  then  reduced  using  the  equal-change-in-velocity 
method  at  five  to  thirty  percent   change  intervals,  providing  approxi- 
mations to  the  original  profile.   Each  of  these  simplified  profiles  was 
used  as  the  environmental  input  to  STUTRACK  I  and  the  source  position 
was  calculated  with  the  transit  times,  entry  angles,  and  horizontal 
ranges  from  the  original  profile. 
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The  position  of  the  source  as  calculated  frora  the  reduced  profiles 
was  compared  to  that  of  the  original  profile.   The  standard  error  of  es- 
timate for  each  of  the  reduced  profiles  was  determined  for  the  36 
profiles  describing  the  twelve  historical  profiles  at  Dabob  Bay.   A 
comparison  of  the  minimum  and  maximum  profiles  of  a  particular  type  with 
the  mid  profile  was  also  performed.   This  comparison  was  designed  to  show 
the  effects  of  using  a  single  profile  to  describe  the  environment  when 
the  actual  profile  existing  at  that  time  is  within  the  range  of  a  parti- 
cular sound  velocity  type.   Table  VIII  shows  the  results  of  these 
comparisons  for  Profile  Type  Four.   All  profile  types  showed  similar 
results . 

In  conducting  the  above  comparisons,  it  was  noticed  that  much  larger 
errors  occurred  at  shallower  depths.   Each  profile  type  shows  its  great- 
est variation  near  the  surface  due  to  water  runoff,  solar  insolation, 
and  seasonal  temperature  variation.   Therefore,  it  was  decided  to  perform 
the  same  analysis  as  above,  but  with  the  data  points  below  the  point 
where  large  variations  in  the  gradients  of  the  minimum  and  maximum  pro- 
files of  a  sound  velocity  type  occurred.   The  calculations  were  performed 
between  100  and  300  feet.   Table  IX  shows  the  results  of  these  calcula- 
tions for  Profile  Type  Four  as  before.   The  values,  as  calculated  from 
data  points  between  100  and  300  feet,  generally,  showed  less  standard 
error  than  the  values  calculated  using  the  same  number  of  points  between 
the  surface  and  200  feet.   Table  X  shows  this  comparison. 
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IX.   CONCLUSIONS 

A.  PRECISION  OF  DATA 

The  analysis  of  the  experiment  showed  conclusively  that  although  the 
magnitude  of  the  spatial  variation  grows  as  the  source's  range  from  the 
array  increases,  the  variation  in  the  angular  description  of  source 
position  is  independent  of  range.   This  is  consistent  with  the  results 
of  the  preliminary  investigation.   It  is  concluded  that  the  increase  in 
the  spatial  variation  of  the  range  position  discrimination  of  the  2.5-MHz 
clock  in  measuring  the  arrival  time  of  the  75-kHz  source  signal  used  in 
the  experiment.   This  uncertainty  can  cause  a  maximum  error  of  3.3  ^j/sec 
in  the  arrival  time  which  is  one-fourth  of  the  period  of  the  75-kHz  sig- 
nal, causing  an  error  of  approximately  one  foot  at  station  1  and  4.5 
feet  at  station  6. 

B.  ENVIRONMENTAL  EFFECTS  ON  RANGE  PRECISION 

The  investigation  of  precision  of  range  position  due  to  the  amount  of 
environmental  information  used  in  the  calculations  has  shown  that  the 
method  of  describing  the  existing  profile  is  not  critical.   Of  the  three 
methods  of  reducing  the  number  of  segments  describing  the  existing  sound 
velocity  profile,  the  equal-change-in-velocity  method  appears  to  produce 
results  with  more  precision  than  either  the  equal-change-in-gradient  or 
"eyeball"  fit  methods.   Describing  a  profile  in  terms  of  a  20  percent 
equal  change  in  velocity  between  segments,  usually  six  or  seven  segments, 
produces  the  best  precision  with  the  minimum  amount  of  environmental 
information  required  for  calculations.   This  result  is  independent  of 
depth  and  horizontal  range. 
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C.   HISTORICAL  ENVIRONMENTAL  ANALYSIS 

The  analysis  of  the  twelve  profile  types  which  describe  the  acoustic 
environment  at  Dabob  Bay  show  that  the  use  of  historical  information 
rather  than  real-time  environmental  data  is  indeed  feasible.   The  use 
of  a  profile  which  has  been  reduced  by  a  20  percent  change  in  velocity 
in  the  equal-change-in-velocity  method  provides  the  most  precision  for 
all  velocity  profile  types  studied.   The  use  of  the  mid  profile  to  de- 
scribe any  profile  falling  within  the  range  of  a  particular  type  causes 
very  little  degradation  in  precision.   Therefore,  it  is  concluded  that 
the  use  of  historical  data  is  feasible  in  the  determination  of  the  source 
position  and  only  a  small  number  of  segments  are  required  to  describe  the 
existing  sound  velocity  profile. 
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X .   COMMENTS  AND  RECOMMENDATIONS 

This  thesis  has  verified  the  preliminary  conclusions  reached  in 
Reference  1,  that  the  precision  of  the  range  position  in  calculations  is 
a  function  of  the  rise-time  of  the  75-kIIz  acoustic  signal  used  in  the 
experiment.   The  angle  of  entry  of  the  signal  to  the  array  is  thus  the 
critical  factor  in  range  precision.   Increasing  the  accuracy  of  time 
measurement  and  decreasing  the  error  due  to  rise-time  (increasing  the 
array  baseline)  will  improve  precision. 

The  historical  profiles  may  be  segmented  as  described  in  this  thesis 
and  catalogued.   The  sound  velocity  profile  which  actually  exists  on  the 
range  could  be  sampled  daily  with  a  velocimeter.   This  profile  can  then 
be  fitted  to  an  historical  profile  type  which  describes  it  and  the  ap- 
propriate catalogued  sound  velocity  profile  used  as  the  environmental 
input  for  any  range  positioning  calculations  to  be  performed  that  day. 

As  noted  in  the  preliminary  investigation,  the  sum  of  the  source 
depth  and  the  individual  calculated  Z  component  at  each  time  differed 
from  the  known  array  depth.  The  range  of  this  difference  was  +4.3  to 
-18.8  feet  with  the  average  being  3.1  feet  less  than  the  array  depth. 
Verification  of  the  recorded  tidal  levels  with  the  predicted,  tides  at 
the  time  the  experiment  was  conducted  showed  no  discrepancies  of  this 
magnitude. 

Further  study  should  be  conducted  to  determine  the  accuracy  of  the 
ranging  system  at  Dabob  Bay  and  Nanoose  and  to  provide  a  greater  con- 
fidence level  in  the  use  of  less  environmental  data  in  range  position 
calculation.:-.   The  apparent  inaccuracy  in  the  Z  component  should 
receive  special  attention. 
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FIGURE    1.       DAB0I3   DAY    RANGE 
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FIGURE    2.       EXPERIMENTAL   AREA 
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FIGURE    5.       RAY    PATH   REFRACTION 
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COMPONENT 
THETA 


STANDARD    ERROR    OF    ESTIMATE 
(DEGREES) 


I                       I                                                DEPTH(FEET)                                                 I 

I                      I         25      T         50       I         75       *       100      I       150       I      200       I 

I    S    '. 
.     T     , 
I    A     1 
!    T 
.     I     . 

."  o  ! 

1    N    ', 

!  1  1  .Cb6  I  .070  I  .060  I  .090  I  .074  I  .078  .* 
!  2  !  .070  I  .036  I  .098  I  I  I  I 
!  3  I  .059  I  .099  I  .112  I  .065  I  .037  !  .102  I 
I  4  I  .131  I  .035  I  .059  J  .053  I  .049  ^  .028  I 
I       5    I       .146    .       .157    .       .038    I       .042    I       .028    "       .037    I 

!       6    I       .092    I       .092    I       .143     .       .078    .       .104    I       .095    ." 

!    5P    I       .124    I       .042    .       .054    !       .051    .       .058    I       .041    I 

',    4P    I       .122    I       .069    I       .042     .'       .026    I       .039     I       .069    .° 

!  3P  I  .045  ."  .030  I  .031  I  .027  I  .061  I  .046  . 
!     IP    !       .055    .*       .046    I       .052    I       .050    I       .066    J       .079    .* 
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COMPONENT 
PHI 


standard  error  of   estimate 
(degrees; 


I                      '.                                               DEPTH(FEET)                                               I 

I                      I         25       I         50       I          75       I       100      '.       150       1       200       I 
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!       4    I       .031    !       .034    !       .037    I       .034    I       .034    \       .035    I 
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!       6    .       .062    I       .055    I       .078    !       .080    ".       .096    .       .080    I 

!  5P  !  .039  I  .037  1  .035  .  .037  *.  .036  *.  .044  I 
!    4P    I       .038    I       .033    I       .037    I       .035    I       .052    I       .032    I 

!  3P  I  .032  '.  .032  I  .030  I  .033  I  .049  I  .043  I 
!     IP    I       .028    ."       .027    I       .033    I       .048    I       .030    '       .056    '. 
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COMPONENT 
R 


STANDARD  ERROR  OF  ESTIMATE 
(FEET) 


'.                       '.                                                DEPTH(FEET)                                                 . 

I                      !         25      .'         50      ! 
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!  a  ! 
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!  n  ! 
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!     1.700    I                    I                    I                    I 
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4     I 
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1       .232    .       .277    ".       .375    1       .338    I 

!    4P    .'       . 
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1       .230    I       .228    I       . 280    I       . 199    I 

!    3P    I      . 

,  165    I       . 195    ! 

!       .178!       .165!       .202!       .166*. 

'   1  p  I     , 

,190    I       .155    , 

!       .167     .       .253    !       .211     .*       -284    I 
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STANDARD  ERROR  OF  ESTIMATE 
CAUSED  dY  ENVIRONMENTAL  DESCRIPTION 


BY  VELOCIMETER 


.  s  , 

.  T  , 

.  A  , 

',      D   1        COMPONENT        '. 

E c 

,     P 

,   T   .  THETA  .   PHI   .    R 

,   H   .  (DEG)  .  (DEG)  .   (FT)  . 

I  1  ! 

!   25  I   .067  I       .0^9  I       .187  I 

50  !   .073  \       .033  \       .114  ! 
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150  "   .076  \       .030  I   .248  I 

!  200  I   .082  ".   .031  \       .175  I 

1  6  ! 

\       25  I   .095  I   .055  '.       .283  ! 

!   50  '.       .095  ™   ,05o  I       .319  I 
!   7r>  '.       .145  I   .079  I   .361  I 

1  100  I   .081  ',       .076  I   .422  I 

!  150  I   .108  I   .098  .'   .415  ! 

!  200  I   .098  I   .081  !   .377  ! 
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STANDARD  ERROR  OF  ESTIMATE  OE  R 
CAUSED  BY  ENVIRONMENTAL  DESCRIPTION 


BY  THE  EQUAL  CHANGE  IN  VELOCITY  METHOD 


.  s  , 

.  T  , 

.  A  , 

!  D  I  PERCENTAGE  CHANGE  ! 
,   E   .                 IN  VELOCITY 

p 

,   T   . 

,   H   .   5    .   10   .   15   .   20   .   25   .   30   . 

I  i  ! 

!   25  !  .153  .  .166  I  .157  '.  .165  .  .135  I  .158  ! 

!  50  I  .142  ."  .151  .'  .149  '.  .129  I  .214  .'  .223  '. 
!   75  '.     .215  1     .140  I  .157  I    .154  .  .142  I    .155  I 

!  100  I  .156  I  .174  1  .165  I  P13d  I  .154  ."  .141  I 

150  I  .195  I  .157  I  .162  ."  .172  '.  .172  !  .158  ', 
!  200  I  .145  I  .lo3  I  .163  l    .130  1  .146  !  .169  '. 

.  6  ! 
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!   75  \    .268  T  .257  I  .314  ',     .580  I  .  i  1  9  .*  .417  I 

\  100  1  .314  *  .246  I  .592  .  .214  I  .834  !  .383  I 
!  150  I  .415  !  .478  I  .274  .*  .314  '.     .319  I  .308  I 

!  200  I  .377  1  .315  I  .319  .  .289  I  .242  1  .984  . 
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STANDARD  ERROR  OF  ESTIMATE  OF  R 
CAUSED  BY  ENVIRONMENTAL  DESCRIPTION 


BY  THE  EQUAL  CHANGE  IN  GRADIENT  METHOD 


.  s 

.  T  ■ 
.  A  , 

,       0       .  PERCENTAGE  CHANGE   . 

.   E   .      IN  GRADIENT 

,   p 

.   T   . 

.   H   .    5   .   10   .   15   . 

I  1  ', 

!   25  I  .170  I  .168  ".  .637  ! 

!  50  1  .139  ".  .151  \  .637  \ 
\       75  \    .146  \     .103  1  .209  \ 

\    100  1  .201  ".  .15  7  ".  .169  \ 

\  150  1  .218  '.  .187  \  .189  1 
\    200  1  .162  ".  .138  1  .177  ! 

'.     6  '. 

25  I  .413  "  .428  I  .435  ! 
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!  150  ."  .310  I  .305  '.     .299  1 

!  200  l    .334  I  .314  I  .305  '. 
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STANDARD    ERROR    OF    ESTIMATE    OF    R 
CAUSED    BY    ENV1RCNMENTAL    DESCRIPTION 


BY    THE     'EYEbALL1    METHOD 


.  s 

.  T  . 

.  A  . 

,        0        , 
.   E   . 

P 
.   T   , 
,   H   , 

!          NUMBER  OF          I 
SEGMENTS 

[    5   ! 

[    3   ! 

!   2  ! 

!   i  I 

I  1  '. 

!   25  ! 

!  .176  ! 

.350  '. 

.  172  ! 

,  .196  '. 

!  50  ! 

.198  ', 

!  .167  ! 

!  .172  ; 

!  .204  I 

',       75  ! 

!  .162  ! 

!  .144  ! 

!  .  130  ', 
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'.   ioo  ! 

!  . 134  ! 

.  .  14o  . 

',     .  132  , 
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!  200  ! 
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\     .140  ! 
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[               NUMBER  OF  SEGMENTS     " 

'.          6   ! 

\          4   1 

1    3 

'.          1   1 

\    6  ', 

!   25 

'.     .241 

'.     .240 

1  .210 

6      *i~  *i*      4 

!  50  ! 

'.     .239  ! 

!  .187  ! 

\     .267 

5£  5£ 

!   75  ! 

!  .118 

'.     .179 

1  .179 

!  ioo  ! 

!  .229  ! 

1  .277  ! 

1  .235  - 

1  .2^5  '. 

!  150 

!  .298  ! 

!  .230  ! 

',     .209 

1  .408  ! 

!  200 

!  .441 

'.     .207 

,   ## 

\     .210  ". 

■*        INDICATES    STANDARD    ERROR 
GREATER    THAN    1.000 
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COMPARISON  OF  STANDARD  ERRORS 
CAUSED  BY  THE 
APPROXIMATION  OF  SOUND  VELOCITY  PROFILES 
BY  THE  EQUAL  CHANGE  IN  VELOCITY  METHOD 


STANDARD  ERROR 

1  1     RANGE  OF  CALCULATIONS 

!  1  0    -  200  FT  \    100  -  300  FT  1 

PER  . 

•  CENT 

!    CHANGE  !    PHI       !  R         !     PHI       *  R 

.(DEG)     .       (FT)     .(DEC)     .       (FT)     . 

*  5  I  0.026*  1.6291  0.031*.  1.684! 
!  10  I  0.031  1  1.98  8 !  0.031*.  1.8441 
!  15  *  0.0161  3.306*.  O.Oll!  1.451* 
!  20  !  0.04l!  5.17q!  0.037!  1.66l! 
!  25  !  0.0421  6.086!  0.0341  1.96l! 
!         30  !    0.023!       6.213!    0«0i2l       2.048! 
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